Spin-Hall transport of heavy holes in III-V semiconductor quantum wells 
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We investigate spin transport of heavy holes in III-V semiconductor quantum wells in the presence 
of spin-orbit coupling of the Rashba type due to structure-inversion asymmetry. Similarly to the case 
of electrons, the longitudinal spin conductivity vanishes, whereas the off-diagonal elements of the 
spin-conductivity tensor are finite giving rise to an intrinsic spin-Hall effect. For a clean system we 
find a closed expression for the spin-Hall conductivity depending on the length scale of the Rashba 
coupling and the hole density. In this limit the spin-Hall conductivity is enhanced compared to its 
value for electron systems, and it vanishes with increasing strength of the impurity scattering. As 
an aside, we also derive explicit expressions for the Fermi momenta and the densities of holes in 
the different dispersion branches as a function of the spin-orbit coupling parameter and the total 
hole density. These results are of relevance for the interpretation of possible Shubnikov-de Haas 
measurements detecting the Rashba spin splitting. 



I. INTRODUCTION 

In the last years, the emerging field of spin electron- 
ics has become a major branch of solid state physics and 
comprises by now all kinds of spin-dependent phenom- 
ena in semiconductor structures and related systems 1-3 . 
Most recently, the possibility of spin-Hall currents has 
attracted considerable theoretical interest 4-19 . In these 
studies a spin current (as opposed to a charge cur- 
rent) driven by an electric field perpendicular to it was 
investigated, where the spinful intinerant charge car- 
rriers are bulk valence-band holes in III-V zinc-blende 
semiconductors 4-9 or conduction band electrons in quan- 
tum wells 10-19 of the same type of materials. In the 
present work we extend these studies to the case of 
heavy holes in quantum wells being subect to spin-orbit 
coupling of the Rashba type due to structure-inversion 
asymmetry 20-24 . 

From a historical perspective, the notion of the spin- 
Hall effect in systems of itinerant spinful charge carri- 
ers was considered first by Dyakonov and Perel 25 in the 
early seventies, and in a more recent paper by Hirsch 26 . 
In these studies the predicted spin-Hall effect is due to 
spin-orbit effects influencing scattering processes upon 
static impurities. Following the terminology used in 10,5 
this is referred to as the extrinsic spin-Hall effect since it 
necessarily requires spin-dependent impurity scattering. 
This is in contrast to the intrinsic spin-Hall effect studied 
theoretically very recently in 4-19 , which is entirely due 
to spin-orbit coupling terms in the single-particle carrier 
Hamiltonian and occurs even in the absence of any scat- 
tering process. We note that this distinction between 
intrinsic and extrinsic effects becomes ambiguous in the 
limit of weak spin-orbit coupling when life time effects of 
carrier quasiparticles have to be taken into account 11 . 

Yet another type of spin-Hall effect was studied re- 
cently by Meier and Loss 27 in a two-dimensional Heisen- 
berg model consisting of non-itinerant spins, in contrast 
to the itinerant-carrier systems mentioned before. 



In the case of conduction-band electrons in III-V semi- 
conductor quantum wells, the intrinsic spin-Hall effect 
results from spin-orbit coupling of the Rashba type (due 
to structure-inversion asymmetry 28 ) and/or the Dressel- 
haus type (due to bulk-inversion asymmetry 29 ). More- 
over, the interplay of these two effects leads to partic- 
ularly interesting transport effects 30 ' 31 ' 12 . For bulk va- 
lence band holes the effects of spin orbit coupling un- 
derlying the intrinsic spin-Hall effect are incorporated 
in Luttinger's effective Hamiltonian 32 giving rise to two 
different dispersion branches differing in their effective 
masses and known as heavy and light holes. Within this 
description, these states form for a given wave vector 
multiplets of total angular momentum 3/2 resulting from 
the s = 1/2 spin degree of freedom of the holes and the 
I = 1 orbital anglar momentum of the p-type atomic or- 
bitals forming the valence band 33 . In a quantum well, a 
splitting between the heavy and light holes occurs due to 
the size quantization in the growth direction of the het- 
erostructure. For a sufficiently narrow quantum well and 
for not too high densities and temperatures, only the low- 
est heavy hole subbands are significantly occupied. Here 
the angular momentum of the heavy hole states points 
predominantly along the growth direction 21 , correspond- 
ing to the total angular momentum states ±3/2, an ap- 
proximation we shall adopt in the present work. 

Experimental investigations of such systems include 
studies of spin polarization and transitions to an insu- 
lating state induced by magnetic fields 34 . Moreover, the 
spin splitting due to spin-orbit coupling has been studied 
in detail via Shubnikov-de Haas oscillations 35 ' 36 including 
also anisotropic properties of the magnetoresistance 37 . In 
recent theoretical investigations the anisotropies in the 
effective ^-factor 38 , spin polarizations 39 , and controlled 
spin rotations induced by Rashba-type spin-orbit cou- 
pling in a spin-FET setup 23 have been considered. In 
the present work we investigate the spin-Hall effect for 
heavy holes in III-V zinc-blende semiconductor quantum 
wells. 
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This paper is organized as follows: In section II we dis- 
cuss the Rashba model for heavy holes in III-V semicon- 
ductor quantum wells. As an aside, we also give explicit 
expressions for the Fermi momenta and the densities of 
holes in the different dispersion branches as a function 
of the spin-orbit coupling parameter and the total hole 
density. Apart from the present investigatios, we expect 
these results to be also of relevance for the interpreta- 
tion of possible Shubnikov-de Haas measurements 40 ' 41 in 
such hole systems. Section III is devoted to the discus- 
sion of our results on spin-Hall transport, and we close 
with conclusions in section IV. 



II. RASHBA SPIN-ORBIT COUPLING FOR 
HEAVY HOLES 

We consider the following single-particle Hamilto- 
nian incorporating spin-orbit coupling due to structure- 
inversion asymmetry for heavy holes in III-V semiconduc- 
tor quantum wells of appropriate growth geometry 20-24 , 
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2m 2Ti 3 V + ; 



(1) 



using the notations p± — p x ±ip y , <r± — a x ±ia y , where p, 
a denote the hole momentum operator and Pauli matri- 
ces, respectively. These Pauli matrices operate on the to- 
tal angular momentum states with spin projection ±3/2 
along the growth direction; in this sense they represent a 
pseudospin degree of freedom rather than a genuine spin 
1/2. In the above equation, m is the heavy-hole mass, 
and a is Rashba spin-orbit coupling coefficient due to 
structure inversion asymmetry across the quantum well 
grown along the [001]-direction chosen as the z-axis. For 
a symmtrically grown quantum well, the coefficient a is 
essentially proportional to an electric field applied across 
the well and therefore experimentally tunable 20-22 . This 
Hamiltonian has two dispersion branches given by 
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where k is the hole lattice momentum and A is the area 
of the system. We note that the validity of above model 
given by the Hamiltonian (1) is restricted to sufficiently 
small wave numbers and densities. In fact, the lower 
of the two dispersion branches (2) is (for a > 0) not 
bounded from below for large wave numbers and has non- 
negative eigenenergies only for k < fi 2 /2ma with a max- 
imum at k = h 2 /3ma. The unbounded decrease of the 
single-particle energies with increasing wave number for 
k > h 2 /3ma is clearly an unphysical feature of the model. 



Therefore, the following considerations are restricted to 
sufficiently small densities such that (at zero tempera- 
ture) only states with wave numbers k < h 2 /3ma are 
occupied. 

For a given Fermi energy Sf and vanishing tempera- 
ture, the above two dispersion branches give rise to two 
different Fermi wave numbers kf fulfilling 
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with k~j < kj (fcj > k f ) for positive (negative) a. Sub- 
tracting these two equations one finds 
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or, cancelling a factor of (fcj~ + kj ), 
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Note that in the relations (5), (6) the Fermi energy does 
not enter explicitly. In fact, the quantity which can be 
immediately controlled experimentally is not the Fermi 
energy but the hole density n given by 
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Combining Eqs. (6), (7), one obtains 
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Note that for fcj < h 2 /3ma the radicand in the above 
equations is always positive. Using 
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one derives from the above equations the following ex- 
plicit expression for k^ as a function of the density n 
and the Rashba parameter a: 
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The difference An of densities of holes in the two disper- 
sion branches, 



(12) 



can be expressed as 
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Such a difference in densities in the two dispersion 
branches could be experimentally detected in Shubnikov- 
de Haas measurements, as it was shown earlier for the 
case of electrons 40,41 . We note that the above results 
for the Fermi momenta and An depend only on the to- 
tal hole densty n and length scale ma/Ti 2 given by the 
Rashba coupling 42 . Winkler et al. 22 have studied both 
theoretically and experimentally the magnitude of the 
Rashba spin orbit coupling in GaAs-based quantum well 
samples with heavy-hole densities of a few 10 14 m~ 2 and 
have found typical values for the characteristic length 
scale ma/h 2 of a few nanometers. In figure 1 we have 
plotted the Fermi wave numbers kf and the difference 
An of holes in the two dispersion branches as a function 
of ma/h 2 at a total hole density of n = 3 • 10 14 m~ 2 . 

On the other hand, solving Eq. (8) for a gives a conve- 
nient expression for the Rashba coefficient as a function 
of the total hole density and the difference of the Fermi 
wave numbers, 
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Finally, writing Eq. (4) in the form 
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one obtains by adding these two equations 
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which does not explicitly depend on the Rashba param- 
eter a. From this relation, it follows with the help of 
Eq. (6) 



Ef = a 
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Using Eqs. (8), (9) it is straightforward to obtain from 
this expression an explicit relation between the Fermi 
energy Ef and the hole density n at a given Rashba pa- 
rameter a. 



III. SPIN-HALL TRANSPORT 

We now investigate the possibilty of spin-Hall trans- 
port of heavy holes in III-V semiconductor quantum wells 
in the presence of Rashba-type spin-orbit coupling. We 
are interested in spin currents (as opposed to charge cur- 
rents) as the linear response of the system to an electric 
field applied in the plane of the well. As we shall see 
below, in such a system the spin current is always per- 
pendicular to the driving electric field and therefore Hall- 
type. We concentrate on the case of zero temperature. 

The single-particle spin current operator is defined by 
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where the factor of 3/2 reflects the angular momentum 
quantum numbers of the heavy holes. The hole velocity 
operator reads v = i\H,r\/K with r being the position 
operator, or, in terms of components, 

'it' h 

A rigorous expression for the spin conductivity, i.e. the 
linear transport coefficient for spin currents driven by a 
spatially homogeneous in-plane electric field, is given by 
the Kubo formula with full frequency dependence for an 
electric field 43 , 



v A(uj + irj) J 

■ J2 f(^(k))(k, M''(t),v y (0)]\k, m) , (21) 

where we have concentrated on the off-diagonal xy- 
componcnts. Moreover, we have assumed zero temper- 
ature T = and non-interacting carriers, which allows 
to formulate the two-body Green's function entering the 
conductivity Kubo formula in terms of single-particle op- 
erators, e = |e| is the elementary charge, and f{e fl {k)) is 
the T = Fermi distribution function for energy e^(k) 
at wave vector k in the dispersion branch a e {+,—}. 
The spin-current operator (in the Dirac picture) for spin 
moment polarized along the ^-direction and flowing in 
the ^-direction is given by 
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jf.»( t ) = eint/H^^-mt/H = ^(t) Px (t) . (22) 

From now on we will assume the Hamiltonian generat- 
ing the above time evolution to include also scattering 
potentials from static random impurities being present 
in the quantum well. The right hand side of Eq. (21) 
is to be understood in the limit of vanishing imaginary 
part r\ > in the frequency argument. This imaginary 
part in the frequency reflects the fact that the external 
electric field is assumed to be switched on adiabatically 
starting from the infinite past of the system, and it also 
ensures causality properties of the retarded Green's func- 
tion occurring in Eq. (21). In general the limiting process 
77 — > does not commute with other limits, and, in par- 
ticular, the dc-limit w — > has to be taken with care 43 . 
In the presence of random impurity scattering, the re- 
tarded two-body Green's function in Eq. (21) will generi- 
cally have a frequency argument with positive imaginary 
part 43 . In this case the limit r\ — > is unproblematic, 
and the imaginary part of the frequency argument is just 
due to impurity scattering and/or other (many-body) ef- 
fects. For the present problem of impurity scattering 
of non-interacting carriers being subject to Rashba-type 
spin-orbit coupling, the resulting imaginary part 77 > 
in the frequency argument is given, to lowest order in the 
Rashba coefficient and the impurity potential, by the in- 
verse of the momentum relaxation time. This is certainly 
a very intuitive result; the formal arguments leading to 
it is completely analogous to the ones used in Ref. 11 and 
can be given along the following lines: In lowest order in 
the spin-orbit coupling and the impurity scattering the 
time-dependent spin-current operator reads 

if ,z (t) « |^S(t)pS(t) (23) 

where the time evolution of (Tq is only due to the Hamil- 
tonian (1), while p x {t) contains the impurity scattering 
but not the spin-orbit coupling. Now it is useful to note 
that, in order to compute the expectation values in the 
Kubo formula Eq. (21), only matrix elements of the time- 
dependent momentum operator p x (t) which are diagonal 
in the wave vector index are needed. This enables to ap- 
ply superoperator techniques developed in Refs. 44 yield- 
ing 

(p£(*))jS « ( e " n °*Px(0))^ - (e- t/T P ° x (0)) m , (24) 

where fio is the scattering master operator in lowest or- 
der of the scattering potential 44 . This operator is the 
same as it occurs as the scattering term in the usual 
Boltzmann equation when evaluated in lowest oder via 
Fermi's golden rule. Thus, Eq. (24) describes the usual 
momentum relaxation due to static impurities in low- 
est order in the scattering potential. For impurity po- 
tentials being isotropic in real space, the momentum p x 
is an exact eigenfunction of f2o, and the eigenvalue is 
given by the well-known inverse momentum relaxation 



time l/r(e) 44 ' 45 which in general depends on the energy 
e(k). To lowest order in the Rashba coupling, this energy 
argument can be replaced with the Fermi energy in the 
absence of spin-orbit interaction. We note that this mo- 
mentum relaxation rate 1/r is the same as obtained in 
the standard diagrammatic approach and thus contains 
the vertex correction for particle transport 43 . However, 
this vertex correction vanishes for short-range isotropic 
scatterers. 

The question of vertex corrections to the spin-Hall 
transport was also examined very recently by Inoue, 
Bauer, and Molenkamp 16 , and by Dimitrova 19 for the 
case of electrons being subject to Rashba spin-orbit in- 
teraction. There the authors reach the conclusion that 
for small but finite disorder the spin-Hall conductivity 
should identically vanish due to vertex corrections, while 
it has its "universal value" of e/87r in the case of a 
perfectly clean system 10-12 . Moreover, Murakami 9 has 
studied vertex correction to spin-Hall transport of holes 
in bulk p-type semiconductors described by a Luttingcr 
Hamiltonian 32 . There the author concludes that vertex 
correction vanish identically, validating, the results of 
Ref. 31 , and ascribes this observation to the fact that the 
underlying Hamiltonian is invariant under inversion of 
momenta, in contrast to the Rashba Hamiltonian. These 
technically rather involved issues are to the opinion of 
the present authors not entirely settled yet and presently 
under investigation. 

Let us now turn to the evaluation of the spin-Hall con- 
ductivity using the aforementioned approximations. A 
straightforward calculation yields 

= -t % -SL f* 7 dk fc4 (25) 

T4mA+ ( w+ i) 2 _ (2a^i) 2 ' 

where, according to the above arguments, the imaginary 
part of the frequency argument is given by the momen- 
tum relaxation rate 1/r. Moreover, the longitudinal spin 
conductivities <J XX Z , a^ z turn out ot be identically zero. 
This is similar to the case of electrons in a quantum well 
being subject to spin-orbit coupling of either the Rashba 
or the Dresselhaus type 10-12 . There a longitudinal spin 
conductivity occurs only if both the Rashba and the Dres- 
selhaus coupling are present 12 . 

The remaning integral in the above expression (25) is 
elementary; however, it leads to a rather cumbersome ex- 
pression which shall not be given here. In the dc limit 
cj = 0, the energy scale of the impurity scattering Ti/t has 
to be compared with the "Rashba energy" er — a(fc") 3 , 

where = y 2mef/h 2 is the Fermi wave number for 
vanishing spin-orbit coupling, which is a typical value for 
k in the inegration in Eq. (25). If the impurity scatter- 
ing dominates over the Rashba coupling, Ti/t » Er, the 
spin-Hall conductivity vanishes with the leading order 
correction given by 
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where the Fermi wave numbers are given by Eq. (11). In 
the opposite case er ^> U/t, the leading contribution to 
the spin-Hall conductivity reads 
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Note that this result for the spin-Hall conductivity de- 
pends only on the length scale ma/fi 2 of the Rashba 
coupling and the total hole density n, but not sepa- 
rately on quantities like the Fermi energy and the effec- 
tive heavy hole mass. If ma/fi 2 is small against the in- 
verse square root of the total hole density (but still fulfill- 
ing er 3> Ti/t), the spin-Hall conductivity approaches a 
value of a^ y z — 9e/87r. This is the case iih/r -C er <C £/• 
This above value should be compared with the univer- 
sal value of e/8n found for the spin-Hall conductivity 
of electrons being subject to a Rashba coupling domi- 
nating possible impurity scattering 10 ' 11 . Thus, in this 
limit, the hole spin-Hall conductivity is enhanced com- 
pared to the case of electrons by a factor of 18, which 
is partially due to the larger angular momentum of the 
heavy holes. In figure 2 we have plotted the spin-Hall 
conductivity for dominating Rashba coupling {er h/r, 
cf. Eq. (27)) as a function of ma/h 2 at a total hole den- 
sity of n = 3 • 10 14 m~ 2 . As seen there, the spin-Hall 
conductivity starts out at cr^y — 9(e/87r) and increases 
with increasing Rashba coupling. 



IV. CONCLUSIONS 

We have studied spin transport of heavy holes in III- 
V semiconductor quantum wells in the presence of spin- 
orbit coupling of the Rashba type due to structure- 
inversion asymmetry. Similarly to the case of electrons, 
the longitudinal spin conductivity vanishes, whereas the 
off-diagonal elements of the spin-conductivity tensor are 
finite giving rise to an intrinsic spin-Hall effect. For a 
clean system we find a closed expression for the spin- 
Hall conductivity depending on the length scale of the 
Rashba coupling and the hole density. In this limit the 
spin-Hall conductivity is enhanced compared to its value 
for electron systems. For dirtier p-doped quantum wells 
when the impurity scattering dominates the spin-orbit 
coupling, the spin-Hall conductivity naturally vanishes 
as also found previously for the case of electrons 11 ' 12 . 
As an aside , we give explicit expressions for the Fermi 
momenta and the densities of holes in the different dis- 
persion branches as a function of the spin-orbit coupling 



parameter and the total hole density. These results are 
expected to be: helpful for the interpretation of possible 
Shubnikov-de Haas experiments aiming at the detection 
of the Rashba spin splitting. 

Note added: After submission of this paper a preprint 
by Wunderlich et aft 6 appeared reporting on an exper- 
imental observation of the spin Hall effect in p-doped 
GaAs quantum wells as studied here. 
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FIG. 1. The Fermi wave numbers k^r and the difference An 
of hole densities in the two dispersion branches as a function 
of the characteristic length ma/h 2 of the Rashba coupling at 
a total hole density of n = 3 • 10 14 m~ 2 . 
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FIG. 2. The spin-Hall condcutivity for dominating Rashba 
coupling {er >• h/r, cf. Eq. (27)) as a function of ma/h 2 at 
a total hole density of n = 3 • 10 14 m~ 2 . 
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